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Abstract. We present a study of instabilities occuring in thick magnetized accretion disks. We calculate the 
growth rates of these instabilities and characterise precisely the contribution of the magneto-rotational and the 
convective mechanism. All our calculations are performed in radially stratified disks in the cylindrical limit. The 
numerical calculations are performed using the appropriate local dispersion equation solver discussed in Blokland 
et al. (2005Jl. A comparison with recent results by Narayan et al. ( 2002J shows excellent agreement with their 
approximate growth rates only if the disks are weakly magnetized. However, for disks close to equipartition, 
the dispersion equation from Narayan et al. 12UU21 loses its validity. Our calculations allow for a quantitative 
determination of the increase of the growth rate due to the magneto-rotational mechanism. We find that the 
increase of the growth rate for long wavelength convective modes caused by this mechanism is almost neglible. On 
the other hand, the growth rate of short wavelength instabilities can be significantly increased by this mechanism, 
reaching values up to 60%. 
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1. Introduction 

Accretion disks are present around a variety of astrophys- 
ical objects, ranging from young proto-stars in star for- 
mation regions to massive black holes in the centers of 
galaxies. A standard model of a geometrically thin accre- 
tion disk was introduced in the early seventies by Shakura 
& Sunyaev (1973). In their model angular momentum is 
transported outwards by an anomalous viscosity mech- 
anism, which is parametrised by the cv-parameter. This 
parameter scales linearly with the radial-azimuthal com- 
ponent of the stress tensor, which is associated with the 
viscous torque providing the angular momentum trans- 
port. The assumption of a geometrically thin disk implies 
that the gravitational energy release by viscous dissipation 
is locally radiated away. 

It was realised in the early nineties that the anomalous 
viscosity mechanism can be provided by the turbulence 
arising from the magneto-rotational instability (Balbus 
& Hawley I1991|) . These authors showed that this essen- 
tially magnetic instability is a very robust one, occuring 
in weakly magnetized thin accretion disks. 

However, in recent years, Chandra observations have 
found examples of undcrluminous black holes at X-ray fre- 
quencies. A good example is our own Galactic Center, Sgr 
A* (Baganoff et al. 12001]) , and the galactic center of the 
elliptic galaxy M87 (Di Matteo et al. ESUOljl . These obser- 



vations might be explained in the context of nonradiat- 
ing accretion flows, which are still using the a description 
from Shakura & Sunyaev 1|1973[) . but are no longer geo- 
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metrically thin (Narayan & LiJTJSEIJ. One of these models 
which has gained a lot of interest in the literature in recent 
years is the convection-dominated accretion flow (CDAF). 
In these type of accretion flow models, the transport of the 
angular momentum outwards by the turbulence arising 
from the magneto-rotational instability is partly counter- 
balanced by transport inwards due to convective motion 
occuring in these thick accretion disks. Therefore the ac- 
cretion rate might be efficiently reduced, which might ex- 
plain the lower observed X-ray luminosities in Sgr A* and 
the central region of M87. 

Recent studies of these CDAF-type of disks, in which 
both convective and magneto-rotational instabilities can 
be found, have been performed by Balbus & Hawley 
( 2002 ) and Narayan et al. ( 2002 ). These authors have used 
a linear stability analysis and identified unstable long- 
wavelength modes with a convective nature, and short- 
wavelength modes with a more dominant MRI behaviour. 

In this paper, we will elaborate on the work from 
Narayan et al. ( 2002 ). These authors have considered a dif- 
ferentially rotating and thermally stratified plasma, with 
a weak axial magnetic field. They use the equation for the 
growth rate as obtained by Balbus & Hawley Q1991[) . who 
used linear stability theory to derive the latter. Narayan et 
al. H2002J) consider five models of an accretion disk, each of 
these five models has a different value of the assumed uni- 



torm Brunt- Vaisala frequency. Iheir results clearly show 
the rise of a plateau in the obtained growth rate at low 
axial wavenumbers as the Brunt- Vaisala frequency is in- 
creased. They argue that the modes in this plateau are 
identified as convective modes, once the Brunt- Vaisala fre- 
quency exceeds the epicyclic frequency, i.e. the H0iland 
criterion (Tassoul. I1987f) . The exact nature of a mode in 
this plateau is identified using results from the linear sta- 
bility theory as performed by Balbus & Hawley 1)19911 
2002). 



We follow up on the above mentioned work, but we 
will numerically solve a sixth order polynomial disper- 
sion relation, using the Local Dispersion Equation Solver 
(LODES), which was already discussed in Blokland et al. 
( 2005 ) . We will not only consider weakly magnetized disks, 
but also consider models which are closer to equipartition. 
We will consider a differentially rotating plasma for which 
the equilibrium quantities are power-law profiles of the 
radius. Furthermore, we explicitly define a value for the 
scale height of the disk H with respect to the radius r, 
with the free parameter e = H/r. 

In our present analysis we only consider axisymmetric 
perturbations, which means that in a model with a purely 
toroidal magnetic field, the magneto-rotational mecha- 
nism is excluded. Therefore in such models, the obtained 
modes are purely convective. The growth rates from these 
modes can be compared to a similar model which includes 
a weak axial magnetic field. This enables one to observe 
the increase of the growth rate due to the presence of the 
magneto-rotational mechanism. 

This paper is organised as follows: in section 2 we dis- 
cuss our model, section 3 shortly recalls the model from 
Narayan et al. (2002), in section 4 we present our results, 
including a comparison with the results from Narayan et 
al. l|2UU2fl . finally in section 5 we present our conclusions. 



2. Accretion disk model 

2.1. Equilibrium of a differentially rotating plasma 

We want to investigate the growth rate of instabilities oc- 
curing in a magnetized accretion disk. In order to quantify 
instabilities using a linear analysis, one has to consider an 
MHD equilibrium. In our case, this equilibrium is one of 
a differentially rotating plasma, which is radially strat- 
ified. In our model the density, pressure, magnetic field 
strength and toroidal velocity only depend on the radius 
r. This type of model is sometimes referred to as an accre- 
tion disk in the cylindrical limit (sec for example Hawley 

We use an equilibrium as in Blokland et al. |2005), 
generalised with an additional parameter a, which is in- 
troduced in order to allow for convective instabilities. The 
following profiles are used for respectively density p, ther- 
mal pressure p, toroidal magnetic field Bg , axial magnetic 



held B z and the toroidal velocity vg: 

p = r-( 3+Q )/ 2 , 
p = e 2 r-( 5+a )/ 2 , 
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in which the a-parameters can be expressed as the ratio 
of both the toroidal and axial magnetic field over the total 
magnetic field B like: 
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and the parameter Vq is defined as: 
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Here, G and M* denote respectively the gravitational con- 
stant and the mass of the central object. Finally, the pa- 
rameter P denotes the radially constant ratio between the 
thermal pressure and the total magnetic pressure B 2 /2, 
i.e. (3 = 2p/B 2 . 

The above power-law profiles satisfy the radial force 
balance equation in cylindrical coordinates: 
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99 , 



(10) 



here g denotes GM,/r 2 , and the prime indicates the 
derivative with respect to the radius r. 

The differences between our equilibrium and the one 
from Narayan et al. (2002) are: 1) we explicitly define the 
above mentioned equilibrium quantities by using power- 
law profiles as a function of radius r; 2) we include the 
toroidal magnetic field component Bg; 3) we define a ratio 
of the scale height H over the radius r, i.e. e = H/r, in 
order to quantify the thickness of our disk model; and 4) 
we do not exclude disks close to equipartition ((3 is a free 
parameter) . 

2.2. Determining the growth rate of instabilities 

We use the local dispersion equation solver (LODES), 
which was recently discussed by Blokland et al. |20O5). 
This code calculates the growth rate of the most unsta- 
ble mode in a given MHD equilibrium for a given radial 
'wavenumber' (q), and a toroidal (m) and axial wavenum- 
ber (k) at a given position r*i. We will only consider ax- 
isymmetric perturbations (m = 0) in this paper. 

In order to determine the radial 'wavenumber' q, we 
use the method as discussed by Blokland et al. (2005). It 



was shown there that under certain assumptions, the nu- 
merical solution of the full set of linearised compressible 
MHD equations governing all MHD modes in disk equi- 
libria obeying Eq. (10), can be avoided for modes obey- 
ing a local dispersion equation found from WKB analy- 
sis. Excellent agreement between the growth rate and the 
eigenfunction behaviour was demonstrated when position 
r\ and associated 'wavenumber' q were properly calculated 
from the full numerical solution. In this paper we follow 
Balbus & Hawley H1991[l and Narayan et al. H2002J) . who 
effectively take q = 0, in order to make a direct com- 
parison with their calculations. We notice that we only 
observed a difference of the order of 1% in our obtained 
growth rates, with respect to our calculations for which 
the 'wavenumber' q was properly calculated. 

The local dispersion equation solver finds the root of a 
sixth order polynomial, using Laguerre's method (Press et 
al. I1988fl . This sixth order polynomial dispersion equation 
represents the local dispersion equation as an approxima- 
tion to the true 10 th degree WKB local dispersion equation 
(Blokland et al. . 120051) . which is more general than the one 
discussed by Narayan et al. ( 2002 ): 1) it allows to calculate 
the growth rate of instabilities for an equilibrium which 
has both an axial and a toroidal magnetic field component; 
and 2) the equilibrium can be either weakly magnetized, 
or close to equipartition. Narayan et al. I|2()()2|) could not 
consider disks close to equipartition, because the fourth 
order polynomial from Balbus & Hawley (1991) was de- 
rived assuming a weak axial and no toroidal magnetic field 
component. 



3. Accretion disks with convection 

3.1. A thick accretion disk 

One of the key parameters of our MHD equilibrium is the 
parameter e, which is taken as a constant free parameter. 
This parameter was introduced by Shakura & Sunyaev 
l|1973f) in their model of geometrically thin accretion disks, 
which are Keplerian rotating. The physical interpretation 
of e in these models is the same as in the MHD equilibrium 
we consider here, i.e. e = c s /ve — H/r, were c s denotes the 
sound speed. MRI instabilities in thin accretion disk mod- 
els were considered by Blokland et al. (|2005|) . However, in 
this paper we consider values of e ~ 1, i.e. we are in the 
regime of sub-Keplerian rotating disks (see e.g. Narayan 
& Yi ll994|) . For these cases we stick to the interpretation 
e ~ H/r, noting that the identification e = c s /vg is not 
valid anymore. Typically, c a /vg in our models will be larger 
up to a maximum factor of order ~ 10. As mentioned in 
section 2, our MHD equilibrium is an example of an accre- 
tion disk model in the cylindrical limit. One has to realise 
that results from these type of models approximate the 
interior of a real height-dependent accretion disk. This 
approximation is correct as long as the axial wavenumber 
k is much larger than the inverse of the scale height H . 
Therefore the identity k » 2n/(er) has to be satisfied in 



1.4 



1 .2 - 



c; 




-1.5 -1.0 -0.5 0.0 

log(kV A /n) 

Fig. 1. Dimensionless growth rate as a function of 
dimensionless wavevector kv^ / fl for accretion disk models 
from Narayan et al. ( 2002 ) . The five growth rates shown 
are from five different solutions: N 2 /Vl 2 ~ 0.0 (solid line); 
N 2 /il 2 = 0.5 (dotted line); N 2 /n 2 = 1.0 (short-dashed 
line); N 2 /{I 2 = 1.5 (long-dashed line); N 2 /tl 2 = 2.0 (dot- 
dashed line). 

our calculations in order to connect these results with the 
interior of an accretion disk. 

3.2. Convective instabilities 

As discussed in section 2, Narayan et al. (2002J) also con- 
sidered an equilibrium of a differentially rotating plasma, 
which only allows for a weak axial magnetic field compo- 
nent. Furthermore, their equilibrium quantities depend on 
radius r, but are not explicitly defined as power-laws obey- 
ing the radial force balance equation (10). Therefore, in 
order to introduce convective instabilities in their model, 
they describe the stratification of the differentially rotat- 
ing plasma in terms of a free parameter N 2 , which is di- 
rectly related to the Brunt- Vaisala frequency Nbv- 

N 2 - -N 2 3 dp dln^/3) 

N ~ NbV ~ hp dr dr ■ (H) 

In a weakly-magnetized rotating medium, the onset of 
convection is determined by the H0iland criterium, i.e. 
TV 2 > k 2 (= 2vs(rvg)' /r 2 ), where n is the epicyclic fre- 
quency. Furthermore, they assume a purely Keplerian sys- 
tem, which means that k 2 = fl^, therefore convective in- 
stabilities will be present in their model when iV 2 > ft 2 . 
Narayan et al. ( 2002 ) present growth rates of instabilities 
as a function of the axial wavenumber. They use the equa- 
tion for the growth rate as derived by Balbus & Hawley 
( 1991J. They normalise their obtained growth rates to the 



rotational frequency, furthermore the wavevector is multi- 
plied by the Alfven velocity Va, and also normalised to the 
rotational frequency. In this way the only free parameter 
left for the growth rate profiles is the parameter N 2 /VI 2 . 
Figure 1 is showing these scaled growth rates for the dif- 
ferent values of N 2 /VI 2 as considered by Narayan et al. 

mm . 

In our cylindrical accretion disk model with power-law 
equilibrium profiles one can calculate N 2 /Vl 2 , and make a 
direct comparison with the work of Narayan et al. (2002). 
Furthermore, it can be shown that in our MHD equilib- 
rium the equality k 2 = £1 2 is always valid, even if the sys- 
tem is not Keplerian rotating (i.e. when e ~ 1). Therefore 
our MHD equilibrium will also be convectively unstable 
once TV 2 > Ct 2 , like in the model from Narayan et al. 
( 2002 ) , the difference being that our models captures sig- 
nificant deviations from Keplerian rotating disks by means 
of the e parameter (fl 2 ^ ^k). 

The H0iland criterium is valid for a hydrodynamical 
equilibrium, but strictly speaking it is not valid for an 
MHD equilibrium which is weakly magnetized. However, it 
can be used as an indication for convective instabilities at 
those wavenumbers where the restoring magnetic tension 
force is much smaller than the buoyancy force (see also 
Christodoulou et al., 2003). These last authors derived a 
stability criterium for axial perturbations for magnetized 
equilibria with purely toroidal magnetic fields. This cri- 
terium is valid for all values of the plasma parameter (3, 
which therefore includes equilibria which are weakly mag- 
netized or close to equipartition. Therefore this criterium 
could be used to determine the presence of convective in- 
stabilities in our models with purely toroidal magnetic 
fields. The criterium can also be derived from the local dis- 
persion equation (see e.g. Keppens et al. (|2002|) ; Blokland 
et al. (E003; Wang et al. QZGGty ). 

3.3. A model for the interior of a thick accretion disk 

In our MHD equilibrium, the parameter e is a free param- 
eter to be interpreted as related to the scale-height H at 
each radial position r. In order to perform a meaningful 
stability calculation for an accretion disk, the wavenum- 
bers of the unstable modes considered must be such that 
they are able to manifest themselves in the interior of the 
accretion disk. These instabilities fit into the interior of the 
considered accretion disk model once the corresponding 
axial wavenumber k > fc m m = 2ir/H. The analysis from 
Narayan et al. i|2002[l shows a maximum of the growth 
rate at log(fcVA/^) ~ 0.0 (see Figure QJ. We therefore 
choose our MHD equilibria such that the numerical value 
of log(fc m i n VA/^) is typically less then ~ —1.5. For these 
models, all unstable modes found with an axial wavenum- 
ber k > & m i n fit into the interior of the considered disk 
model, and include the most unstable ones. 

We found that in order to satisfy the restriction k > 
fcmin with log(fc m ;„VA/r2) ~ —1-5 a high ratio of the 
toroidal magnetic field strength to the axial magnetic 
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Fig. 2. Logarithmic gray-scale plot of the value 
log(fc m ; n VA/^) in the model for which a± = 300.0, c%2 = 
1.0 and a = —3.0. The solid thick line corresponds with 
N 2 — il 2 . The thick solid curve for the H0iland criterium 
deviates from a straight vertical line, because of the dy- 
namical importance of the magnetic field on the equilib- 
rium for low values of (3. 



field strength is needed. Figure [21 shows the parameter 
log ( fc m i n Va/Q) as a function of the parameters e and (3 
of our model, for which the other parameters have been 
fixed. Indeed for most values of e and (3 the correspond- 
ing numerical value log (fc m in Va/Q) < —1.5. The thick 
solid line in Figure marks the H0iland criterium, i.e. 
N 2 = fl 2 . The models on the right hand side of this line 
are convectively unstable, the models on the left hand side 
are convectively stable. Notice that the line N 2 = fi 2 can 
be drawn uniquely as a function of (3 and e once ot\ and 
«2 are determined. This is because the numerical value 
of N 2 /H. 2 , calculated from the equilibrium presented in 
section 2, does not depend on the radius r. 

As argued above we will choose a configuration of an 
accretion disk with most of the magnetic field strength in 
the toroidal component, in order to allow for both con- 
vective and magneto-rotational instabilities. The values 
of the ratio of the Brunt- Vaisala, frequency and the rota- 
tional frequency marks the free parameter in the paper 
from Narayan et al. ( 2002 ) . We plot the numerical value of 
this parameter in our models in figure El We have plotted 
contours with N 2 /VI 2 = 0.5, 1.0, 1.5, 2.0, corresponding to 
the models discussed in Narayan et al.( 2002). It can be 
noted that in very weakly magnetized, thick disk models 
we obtain high values for N 2 /ft 2 , the highest value being 
N 2 /fl 2 = 600 (model A in Table 1). Those are expected to 
be rigorously convective in a pure hydrodynamical sense. 




0.7 0.8 0.9 1.0 



Fig. 3. Gray-scale plot of the parameter 7V 2 /f2 2 in the 
model for which a-y = 300.0, a 2 = 1-0 and a — —3.0. 

4. Accretion disks with two instabilities 

4.1. Introduction 

We present the results from a linear analysis performed 
for different MHD equilibria, using the local dispersion 
equation solver (LODES). In our calculations we will only 
consider equilibria with a constant density (a = —3.0), 
while most of the magnetic held strength is in the toroidal 
direction (a\ = 300.0; a 2 = 1.0). Finally, we will only 
consider thick accretion disk models by considering MHD 
equilibria with e > 0.8, in order to allow for mixed convec- 
tive magneto-rotational instabilities. We will consider five 
different MHD equilibria and calculate the growth rates 
of the unstable modes. Additionally, we will repeat these 
calculations for all five models with the only difference 
that the axial magnetic held equals zero (a 2 — 0.0). Since 
we are only considering axisymmetric perturbations, these 
last calculations yield instabilities which can not have a 
magneto-rotational nature, and are purely convective in- 
stabilities. Because the axial magnetic helds are weak, 
comparing the growth rates of a model with or without 
the axial magnetic held component provides a quantita- 
tive measure of the magneto-rotational contribution of an 
unstable mode in the interior of an accretion disk. 

4.2. Weakly magnetized disks: from sub-Keplerian to 
near-Keplerian disks 

In this subsection, we will investigate the influence of the 
parameter e on the nature of the instabilities in weakly 
magnetized {j3 = 1000) MHD equilibria (models A-D, 
Table 1). We will present the growth rates of the instabil- 
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0.88 


1000.0 
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0.845 
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12.0 
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1.0 


10.0 
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Table 1. Parameters of accretion disk models: all the 
models considered have oti = 300, a = —3 and = 1 
(a 2 = 0) for models with (without) axial magnetic field. 



ities in our model as a function of the axial wavenumber, 
using the same scaling as Narayan et al. ( 2002 ) . 

Figure 0] shows our results for model A, which has a 
corresponding value of the ratio N 2 /fl 2 = 600 (model 
A, Table 1). The left panel shows a comparison between 
our calculations for the growth rate of the instabilities 
(crosses) and the calculation from Narayan et al. (2002) 
(solid line). There is perfect agreement between both re- 
sults. The right panel shows the ratio of the growth rate 
of model A with an axial magnetic field (a 2 = 1.0) 
with respect to model A without an axial magnetic field 
(a 2 = 0.0). The model without axial magnetic field ex- 
cludes the magneto-rotational mechanism completely. One 
can see that the ratio never exceeds unity, implying that 
there is no enhancement of the growth rate due to the 
magneto-rotational mechanism. Furthermore, we observe 
in the left and right panel that there is no observed max- 
imum value of the growth rate at log(£;i>A/^) w 0.0, a 
feature also associated with the magneto-rotational insta- 
bility. Therefore, the observed instabilities have a (nearly) 
completely convective nature in this particular disk model, 
as expected from the very high value of N 2 /£l 2 . 

Figure shows our results for model B, which has a 
corresponding value of the ratio N 2 /£l 2 — 4.5 (model B, 
Table 1). One can see in the left panel that the growth rate 
from our calculations (crosses) again matches the solution 
for the growth rate (solid line) from Narayan et al. ( 2002 ) , 
despite the fact that we consider an overall weak, but pre- 
dominantly toroidal magnetic field. The right panel of fig- 
ure El shows the ratio of the growth rate of model B with 
axial magnetic field (a 2 — 1.0), with respect to model B 
without an axial magnetic field component (a 2 = 0.0). 
It is observed that there is an overall small contribution 
(less than 1 percent) to the growth rates by the magneto- 
rotational mechanism. This indicates that the unstable 
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Fig. 4. Left panel: dimensionless growth rate as a function of dimcnsionless wavenumber Uva/^ for MHD equilib- 
rium A (see table 1). The solid line corresponds with the solution from Narayan et al. I|1994f) . The crosses correspond 
to the solution using the code LODES. Right panel: ratio of the growth rate 7mc with (a.2 = 1-0) and the growth rate 
7c without {a.2 = 0.0) an axial magnetic field component as a function of the dimensionless wavenumber kvp^/Sl. 
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Fig. 5. Left panel: dimensionless growth rate 7/fi as a function of dimensionless wavenumber kvj^/Vt for MHD equilib- 
rium B (see table 1). The solid line corresponds with the solution from Narayan et al. I|1994f) . The crosses correspond 
to the solution using the code LODES. Right panel: ratio of the growth rate 7mc with (a.2 = 1.0) and the growth rate 
7c without (o!2 = 0.0) an axial magnetic field component as a function of the dimensionless wavenumber kvx/Sl. 



modes have a dominant convective nature in this model, 
like in model A. 

In model C, the thickness of the disk model has been 
reduced further, such that the obtained ratio N 2 /Q 2 ~ 2.0 
equals the maximum value Narayan et al. I|2002|l have con- 
sidered in their paper. The left panel of Figurc|Slshows the 
calculated growth rates (crosses) of this MHD equilibrium, 
and there is a perfect match with the solution (solid line) 



from Narayan et al.|2002). The right panel of Figure El 
shows the ratio of the growth rate of model C with axial 
magnetic field (02 = 1.0) with respect to model C without 
an axial magnetic field component (012 = 0.0). The contri- 
bution from the magneto-rotational mechanism is signifi- 
cant, of order 25 percent around values of the wavenumber 
where the growth rate has its maximum: this increase in- 
dicates that the unstable modes in this range of wavenum- 
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Fig. 6. Left panel: dimcnsionless growth rate as a function of dimcnsionless wavenumber kvx/Q for MHD equilib- 
rium C (see table 1). The solid line corresponds with the solution from Narayan et al. I|1994f) . The crosses correspond 
to the solution using the code LODES. Right panel: ratio of the growth rate 7mc with (a.2 = 1.0) and the growth rate 
7c without {a.2 = 0.0) an axial magnetic field component as a function of the dimensionless wavenumber kvp^/Sl. 
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Fig. 7. Left panel: dimensionless growth rate 7/fi as a function of dimensionless wavenumber kvj^/Vt for MHD equilib- 
rium D (see table 1). The solid line corresponds with the solution from Narayan et al. I| 1994(1 . The crosses correspond 
to the solution using the code LODES. Right panel: ratio of the growth rate 7mc with (a.2 = 1.0) and the growth rate 
7c without (a.2 = 0.0) an axial magnetic field component as a function of the dimensionless wavenumber kv\/Sl. 




bers are significantly amplified by the magneto-rotational 
mechanism. Furthermore, model C shows that for axial 
wavenumbers \og(kvx/VL) < —1.0, the contribution of the 
magneto-rotational mechanism has decreased to a few per- 
cent. This indicates that the unstable modes in this range 
of wavenumbers have a dominant convective nature. This 
conclusion on the long wavelength regime was also made 
by Narayan et al. (2002) in their analysis. 



Finally, Figure [7| shows our results for model D, which 
has a corresponding value of the ratio N 2 /Vl 2 — 1.5 (see 
Table 1). The left panel of Figure [3 shows the calcu- 
lated growth rates (crosses) of this MHD equilibrium, in 
agreement with the solution (solid line) from Narayan 
et al. l|2UU2|l . The right panel of Figure [7| shows the ra- 
tio of the growth rate of model D with axial magnetic 
field (pt2 = 1.0) with respect to model D without an ax- 
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Fig. 8. Left panel: dimensionless growth rate as a function of dimcnsionless wavenumber kvA/^l for MHD equilib- 
rium E (see table 1). The solid line corresponds with the solution from Narayan et al. 1)1994)1 . The crosses correspond 
to the solution using the code LODES. Right panel: ratio of the growth rate 7mc with {a.2 = 1-0) and the growth rate 
7c without {a.2 = 0.0) an axial magnetic field component as a function of the dimensionless wavenumber kvp^/Sl. 
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Fig. 9. Left panel: dimensionless growth rate 7/fi as a function of dimensionless wavenumber kvj^/VL for MHD equilib- 
rium F (see table 1). The solid line corresponds with the solution from Narayan et al. 1)1994(1 . The crosses correspond 
to the solution using the code LODES. Right panel: ratio of the growth rate 7mc with (a.2 = 1.0) and the growth rate 
7c without («2 = 0.0) an axial magnetic field component as a function of the dimensionless wavenumber kv\/£l. 



ial magnetic field component {0L2 = 0.0). The contribu- 
tion from the magneto-rotational mechanism is increasing 
up to 60 percent around values of the wavenumber where 
the growth rate has its maximum. For axial wavenum- 
bers \og(kvA/Q) < —1.0, the contribution of the magneto- 
rotational mechanism has a maximum of order 5 percent, 
which indicates their dominant convective nature. 



The models C and D discussed above have the same 
ratio of iV 2 /fi 2 as the models discussed by Narayan et 
al.(|2D02). Our conclusions confirm those by Narayan et 
al. 1)2002(1 for the unstable long- wavelength modes (&Va <C 
57) in these two models: they are clearly convective modes. 
However, the wavelength modes in these models with 
kVx ~ fi are shown to be maximally amplified by the 
magneto-rotational mechanism: the amplification is in- 



creasing as the ratio TV /Si is decreasing. However, these 
modes seem to be determined by both the convective 
and the magneto-rotational mechanism. Finally, for modes 
with kVA > Si, the amplification of the growth rate is 
shown to be decreasing, due to the restoring magnetic ten- 
sion force. 

4.3. From weakly magnetized disks to disks close to 
equipartition 

In this subsection we investigate the influence of the mag- 
netization on the nature of the instabilities in our consid- 
ered MHD equilibria (models E & F, Table 1). 

Figure|Hlshows our results for model E. The considered 
MHD equilibrium has the same parameters as model A 
from the previous subsection, except this one is closer to 
equipartition ((3 = 20). Due to the stronger magnetic field, 
this equilibrium has a lower value of the ratio N 2 /Si 2 (see 
Table 1). Our calculations (crosses) and the calculation 
from Narayan et al. ( 2002) (solid line) deviate from each 
other, since the equation used by Narayan et al. {2002) 
is not valid for this model. This is because the toroidal 
magnetic field strength can not be neglected anymore in 
the linear analysis, a component which was not taken into 
account in the equation used by Narayan et al. (|2005). 
Contrary to the model A, this model shows a maximum 
value of the growth rate, close to the typical value of the 
magneto-rotational instability, i.e. \og(kvp^/Vl) ~ 0.0. 
The right panel shows the ratio of the growth rate of model 
E with axial magnetic field {012 = 1-0) with respect to 
model E without an axial magnetic field component («2 = 
0.0). There is a slight enhancement observed close to the 
maximum value of the growth rate: this indicates that 
the magneto-rotational mechanism starts to contribute to 
the growth rate of unstable modes in this region. Notice 
that the enhancement is still very weak, of order 1 percent 
(right panel Figure 0. 

Finally, figure U3 shows our results for model F. In this 
MHD equilibrium the magnetization is close to equiparti- 
tion ((3 = 10.0), with a corresponding value N 2 /£l 2 = 6.0 
(see Table 1). We notice that we do not consider a disk 
in exact equipartition (i.e. (3 — 1.0), because this model 
would not allow for convective instabilities (i.e. N 2 /Si 2 < 
1 as can be seen in Figure 3). The left panel of Figure |5] 
shows that the approximation from Narayan et al. (2002) 
(solid line) is no longer valid. Taking into account the 
toroidal magnetic field component, our calculations show 
a distinct maximum which occurs at the typical value 
log(feuA/Sl) = 0.0 for magneto-rotational instabilities. The 
right panel of figure El shows the ratio of the growth rate 
of model F with axial magnetic field [0.2 = 1.0) with re- 
spect to model F without an axial magnetic field compo- 
nent («2 = 0.0). It can be observed that in the region 
around the maximum of the growth rate, the contribu- 
tion of the magneto-rotational mechanism is of order 30 
percent. This increase indicates that the unstable modes 
in this range of wavenumbers are significantly amplified 



by the magneto-rotational mechanism. However, tor axial 
wavenumbers log(fcvA/Sl) < —1.0, the contribution of the 
magneto-rotational mechanism has decreased to a few per- 
cent. This indicates that the unstable modes in this range 
of wavenumbers have a dominant convective nature. 

5. Conclusions 

We calculated growth rates of instabilities for MHD equi- 
libria of magnetized accretion disks. The disks models 
were taken in the cylindrical limit and included both 
toroidal and axial magnetic field components. Due to 
the presence of a toroidal magnetic field component, the 
unstable modes we find are strictly speaking overstable 
modes (see e.g. Blokland et al.. l2005|) . The parameters of 
the model in our calculations were taken such that the 
axial wavenumber of the most unstable modes were larger 
than the inversed scale height H . Therefore the analysis 
can be used as an approximation for a stability analysis 
of the interior of an accretion disk. Our models consid- 
ered both weakly magnetized disks, as well as disks which 
are close to equipartition. All calculations were performed 
with a restriction to axisymmetric perturbations (to = 0). 
Our most important results are summarised below: 

— We have considered models of thick accretion disks, 
which are subject to convective magneto-rotational in- 
stabilities, and we quantified the contribution of the 
MRI mechanism of all the instabilities considered. 

— Our calculations are diverging from the solutions from 
Narayan et al.(2002), when the disk models are be- 
coming significantly magnetized. The deviations are 
contributed to the toroidal magnetic field component, 
which was neglected in the derivation of the equation 
for the growth rate used by Narayan et al.(2002). 

— Our calculations show that the contribution from the 
magneto-rotational mechanism to the growth rate be- 
comes significant as the disks get close to equipartition. 
The contribution peaks for modes near the maximum 
value of the growth rate. 

— Our calculations for the growth rates agree with the 
calculations from Narayan et al. (2002) for weakly 
magnetized accretion disk models, even though we in- 
clude dominant toroidal field components. 

— Our calculations show that the contribution from the 
magneto-rotational mechanism to the growth rate be- 
comes significant as the scale-height H is decreased. 
This contribution peaks for modes near the maximum 
value of the growth rate. 

— All our calculations show that the contribution from 
the magneto-rotational mechanism to the growth 
rate becomes neglible for unstable modes with 
log(/ci>A/Sl) < —1.0, which are safely interpreted as 
purely convective modes in all cases. 
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